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Abstract. We have computed accurate 1-D solar models including both a macroscopic mixing process in the solar 
tachocline as well as up-to-date microscopic physical ingredients. Using sound speed and density profiles inferred 
through primary inversion of the solar oscillation frequencies coupled with the equation of thermal equilibrium, 
we have extracted the temperature and hydrogen abundance profiles. These inferred quantities place strong 
constraints on our theoretical models in terms of the extent and strength of our macroscopic mixing, on the 
photospheric heavy elements abundance, on the nuclear reaction rates such as Su and S34 and on the efficiency of 
the microscopic diffusion. We find a good overall agreement between the seismic Sun and our models if we introduce 
a macroscopic mixing in the tachocline and allow for variation within their uncertainties of the main physical 
ingredients. From our study we deduce that the solar hydrogen abundance at the solar age is Xi nv = 0.732 ±0.001 
and that based on the 9 Be photospheric depletion, the maximum extent of mixing in the tachocline is 5% of the 
solar radius. The nuclear reaction rate for the fundamental pp reaction is found to be Sii(O) = 4.06 ± 0.07 10 -25 
MeV barns, i.e., 1.5% higher than the present theoretical determination. The predicted solar neutrino fluxes are 
discussed in the light of the new SNO/SuperKamiokande results. 
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1. Introduction 

Over the past decade our understanding of the solar in- 
terior has improved significantly. Today with the pre- 
cise helioseismic data available from the GONG (Global 
Oscillation Network Group) ground based instruments 
and the SOHO (SOlar and Heliosphcric Observatory) 
space experiments (Gough et al. 1996; Thompson et 
al. 1996; Frohlich et al. 1997; Gabriel et al. 1997; Schou 
et al. 1998), the detailed internal structure and complex 
dynamics of our star can be inferred with reasonable ac- 
curacy using inversion techniques. In addition to sound 
speed c and density p profiles, the internal rotation rate 
SI can also be inferred. It reveals that on the top of an 
almost uniformly rotating radiation zone (with a rotation 
period of about 28 days), the bulk of the convection zone is 
differentially rotating with properties close to what is de- 
duced from sunspot tracking, i.e., a period at the equator 
of 25 days and at the pole of 33 days corresponding to a 
contrast Af2 of 30%. The sharp transition region between 
these two distinct zones, located around 0.7 Rq, has been 
called the tachocline (Spiegel & Zahn 1992); it is thought 
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to play an important role in determining the structure and 
the chemical evolution of the Sun (Brun, Turck-Chieze & 
Zahn 1999; Elliott & Gough 1999). 

Indeed, it appears that the solar structure deduced 
from helioseismology and the observed photospheric com- 
positions can not be explained adequately without in- 
voking some mixing in the radiative interior (Brun et 
al. 1999). This conclusion has been drawn after a careful 
study of the microscopic processes present in solar mod- 
els (Christensen-Dalsgaard et al. 1996; Morel, Provost & 
Bcrthomicu 1997; Bahcall, Basu & Pinsonneault 1998a; 
Brun, Turck-Chieze & Morel 1998). More precisely, in 
the early 90's, after significant improvements in the de- 
scription of the solar plasma through better equation of 
state, opacities and nuclear reaction rates, helioseismic 
studies have established the need for microscopic diffu- 
sion of helium and heavy elements in the radiative inte- 
rior (Christensen-Dalsgaard et al. 1993). But it was soon 
realized that models including only microscopic diffusion 
exhibit sharp composition gradients below the base of the 
convection zone which are not consistent with helioseismic 
data; these favour instead smoother composition profiles 
within this region (Basu & Antia 1994). In spite of fur- 
ther improvements in solar models, this discrepancy still 
persists around 0.7 Rq (Christensen-Dalsgaard et al. 1996; 
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Bran et al. 1998), suggesting that some extra mixing must 
be implemented in the models. 

Further, using primary inversions for sound speed and 
density and the equations of thermal equilibrium, Antia & 
Chitre (1998) have inferred the hydrogen abundance pro- 
file in the radiative interior. This profile confirms the pres- 
ence of such mixing in the Sun, as the hydrogen abundance 
appears to be almost constant in the region r > O.68i?0. 

Another evidence for mixing occurring in that region 
comes from the photospheric light elements composition. 
Purely microscopic processes cannot reproduce the under 
abundance of lithium observed in the Sun and in open 
clusters (Grevesse, Noels & Sauval 1996; Cayrel 1998; 
Richard et al. 1996; Turcottc et al. 1998; Brun et al. 1999). 

We are thus compelled to introduce some mixing pro- 
cesses in the stably stratified radiative interior. The pos- 
sible causes of instabilities leading to such mixing are the 
solar rotation, the magnetic field or penetrative convec- 
tion (Zahn 1998). The recent study of Balachandran & 
Bell (1998) on the photospheric light elements abundance 
of 7 Li and 9 Be puts strong constraints on the extent, am- 
plitude and location of such instabilities. It is now be- 
lieved that only the lithium is significantly depleted, by 
more than a factor of 100 in comparison to the meteoritic 
composition, while the beryllium has varied by only 10% 
over the last 4.6 Gyr. The temperatures at which these 
two species are destroyed by nuclear burning are respec- 
tively ~ 2.7 x 10 6 K (at O.66i? ) and - 3.2 x 10 6 K (at 
O.59i?0), which are relatively close to the temperature at 
the base of the convection zone ~ 2.2 x 10 6 K at 0.713i?o 
(Christensen-Dalsgaard et al. 1991). This implies that any 
macroscopic processes for such lithium destruction have to 
be located near the top of the radiation zone and cannot 
extend deeper than (~ 8%) in solar radius without pro- 
ducing an excessive destruction of 9 Bc. 

This requirement is satisfied if the mixing is confined 
in the tachocline. For this reason, Brun et al. (1999) 
calculated the mixing occurring in that layer, based 
on Spiegel & Zahn's hydrodynamical description of the 
tachocline. They found indeed that such mixing improved 
the agreement between the models and the Sun, pro- 
vided the secular variation of the tachocline was taken 
into account. Alternative approaches based on gravity 
waves (Montalban & Schatzman 1996) and magnetic field 
(Barnes, Charbonneau & MacGregor 1999) have also been 
studied. 

In this paper we intend to go further in understand- 
ing the influence of tachocline mixing on the solar inter- 
nal structure and composition by using the powerful seis- 
mic diagnostic. In particular, we examine the cumulative 
effect of varying, within their intrinsic uncertainties, the 
basic microscopic ingredients such as the nuclear reaction 
cross sections, the microscopic diffusion coefficients, etc. 
Seismic probing of the solar structure has been improved 
to a level where it can be used to constrain physical pro- 
cesses, since the Sun provides a ready-made cosmic lab- 
oratory for testing various aspects of physics. For exam- 
ple, there have been some attempts to constrain the nu- 



clear reaction rate for pp reaction using hclioseismic data 
(Antia & Chitre 1998, 1999; DegPInnocenti, Fiorentini & 
Ricci 1998; Schlattl, Bonanno k Paterno 1999). These 
studies indicate that the cross-section for pp reaction needs 
to be increased by a few percent over the currently ac- 
cepted value (Adelberger et al. 1998). Using a similar ap- 
proach Weiss, Flaskamp & Tsytovich (2001) found that 
enhancing the electron screening by about 5% improves 
the agreement between solar model and helioseismically 
inferred sound speed. Therefore, we would like to revise, 
among other quantities, the pp cross-section deduced from 
hclioscismology, by using in our solar models cither the 
weak (Salpeter 1954) or intermediate (Mitler 1977) treat- 
ment for electron screening. We also deduce the photo- 
spheric hydrogen abundance AT p h, the maximum extent h 
of the tachocline mixing allowed in a solar model and pre- 
dict the theoretical neutrino fluxes in light of the recent 
SNO results (Ahmad et al. 2001, 2002). 

We have organised our paper as follows. In section 2, 
we briefly recall how we compute our 1-D solar models 
with or without the presence of tachocline mixing, de- 
scribe our inversion techniques of the solar acoustic fre- 
quencies and demonstrate the need for further progress 
in solar modelling. In section 3, we present our latest re- 
sults on the sound speed, density, hydrogen abundance 
and temperature profiles obtained with our modified so- 
lar tachocline models and discuss the resulting neutrino 
fluxes. Finally in section 4, we comment on our findings 
and outline our conclusions. 

2. Modelling Approach 

2.1. Construction of solar models with tachocline 
mixing 

In order to model the Sun, we use the CESAM code (Morel 
1997), which solves the structure equations (Kippenhahn 
& Weigert 1994) in time and space for a spherically sym- 
metric star of one solar mass (Mq) in mechanical and 
thermal equilibrium. Once the evolved structure of a given 
model reaches the age of 4.6 Gyr, including a pre main 
sequence (PMS) phase of ~ 50 Myr, it is calibrated to 
the solar radius Rq , luminosity Lq and surface heavy ele- 
ments abundance (Z/X) p \ 1 to within an accuracy of 10 -5 
(see Table |l|) . This is done by modifying the mixing length 
parameter a, the initial helium Yq and heavy elements 
Zq abundances. This accurate calibration of solar values 
is crucial and allows us to test different solutions in our 
search for the best agreement between our models and the 
Sun. We refer to Brun et al. (1998) for a more detailed de- 
scription of our solar models. 

The macroscopic mixing present in the tachocline is 
modelled by adding to the equation for chemical evolution 
an effective time dependent diffusion coefficient Dr(r, f), 
based on the hydrodynamical description of the tachocline 
developed by Spiegel & Zahn (1992). In their study they 
invoked the anisotropy of the turbulence in a stratified 
medium to explain the thinness of this layer. We refer to 
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Table 1. Solar Observations: physical parameters, helio- 
seismic observations, solar neutrino detections 



Physical parameters 

M© = (1.9891 ± 0.00 04) x 10 33 g 

R e = (6.9599 ± 0.0002) x 10 10 cm 

L Q = (3.846 ± 0.004) x 10 33 erg s" 1 

Age = 4.6 ± 0.04 Gyr 

(Z/X) ph = 0.0245 x (1 =b 0.1) 

Helioseismic observations 

F surf = 0.249 ± 0.003 

Rbcv/Re = 0.713 ± 0.003 

h/Ro < 0.05 (tachocline thickness) 

Solar neutrino detections 

71 Ga = 75 ± 5 SNU (average of all Gallium experiments) 
37 CI = 2.56 ± 0.23 SNU (Homestake) 
H 2 = 2.32 ± 0.08 x 10 6 cm -2 s" 1 (SuperKamiokande) 
D 2 = 1.75 ± 0.14 x 10 6 cm" 2 s" 1 (SNO-charged current) 
D 2 = 5.09 ± 0.62 x 10 6 cm" 2 s -1 (SNO-neutral current) 



Note: For the gallium and chlorine detections of the solar neu- 
trino flux we adopt the standard unit, 1 SNU = 10~ 36 cap- 
tures/atom/s. 

Brun et al. (1999) for a complete description of the differ- 
ent steps followed to deduce from their model the effective 
turbulent diffusivity Dt(t, t) used in this work. This coef- 
ficient depends on two parameters: 

• the tachocline thickness at the solar age, h (or the closely 
related quantity d ~ 2/i), a relatively well-known quantity 
(Antia et al. 1998; Corbard et al. 1999, see Table 1), 

• the Brunt- Vaisala, frequency M/2tt, which varies with 
depth and is taken as constant in this model, representing 
some average over the tachocline; furthermore, it depends 
on the extent of overshoot and its value is therefore some- 
what uncertain. 

The time dependence of the angular velocity is based on 
the Skumanich law (Skumanich 1972), e.g., Q(t) oc t -1 / 2 . 
This law is not adequate for the early phases of the so- 
lar evolution when the star contracts and/or exchanges 
angular momentum with its accretion disk (see Piau & 
Turck-Chieze 2001). Nevertheless, for this study concerned 
mainly with the present day Sun, it is satisfactory. 

2.2. Inversion techniques 

To test and constrain the solar models, we compare their 
sound speed, density, temperature and hydrogen abun- 
dance profiles with seismically deduced ones. The sound 
speed and density profiles are inferred using a Regularised 
Least Squares (RLS) inversion technique (Antia 1996). 
This primary inversion is based only on the equations of 
mechanical equilibrium, and has been tested through ex- 
tensive comparisons (e.g., Gough et al. 1996). For these 
primary inversions, we use a set of modes in the range 
of harmonic degree £ < 190, obtained from the first 360 
days of operation of the Michelson Doppler Imager (MDI) 
(Schou et al. 1998). The inversion results in the central 



region depend on the set of low degree acoustic modes 
used (i.e., £ < 2, Basu et al. 2000) but not to the extent of 
modifying the conclusions of this work. To infer secondary 
quantities such as temperature and chemical composition 
within the Sun, we follow the treatment given by Antia & 
Chitre (1998). 

Apart from evolutionary solar models, we also con- 
struct some static ones, using a composition profile calcu- 
lated by the evolutionary stellar structure code CESAM. 
These models use the same physical inputs as the evolu- 
tionary models, but include a different treatment of at- 
mosphere, using the atmospheric model of Vernazza et 
al. (1981) as well as the opacity tables from Kurucz (1991) 
at low temperatures and the formulation of Canuto & 
Mazzitclli (1991) to calculate the convective flux. Because 
of these differences the surface layers in the Sun are better 
represented in these static models. 

2.3. Earlier results 

Before introducing any new modifications in our evolu- 
tionary solar models, let us recall what are the strengths 
and weaknesses of the tachocline models of Brun et 
al. (1999). 
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Fig. 1. The relative difference in sound speed and density 
profiles between the Sun and solar models is shown as a 
function of fractional radius. Model ref including only mi- 
croscopic diffusion is represented with a solid line and the 
mixed models B tz and its static equivalent with respec- 
tively a short-dashed and long-dashed lines. Superimposed 
on model ref is the la error envelope coming from the 
helioseismic inversion. 

In Fig. |l| we represent the relative sound speed and 
density differences Sc/c and dp/ p between the seismic Sun 
obtained using the inversion procedure described above 
and our purely microscopic diffusive model, hereafter re- 
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ferred to as ref model or a typical tachocline model of Brun 
et al. (1999), namely, model B tz . A first quick look reveals 
that the mixed model B tz shows a better overall agree- 
ment both for density and sound speed compared to the 
reference one. For both these quantities, the transition at 
the base of the homogeneous convection zone is smoother, 
resulting in an almost disappearance of the pronounced 
peak seen in Sc/c of model ref. Further, the hydrogen 
abundance in the convective envelope is slightly closer to 
the seismic one, due mainly to the limiting action of the 
macroscopic mixing on the gravitational settling of the 
chemical elements (see Table 3). Brun et al. (1999) found 
that macroscopic mixing at the base of the convection 
zone reduces by 25% the microscopic diffusion in compar- 
ison to a purely mi cros copic model such as ref. Moreover, 
as we shall see in §3^3, the hydrogen abundance profile is 
smoother and does not exhibit any sharp gradient just be- 
low the base of the convection zone. Further, the lithium 
depletion achieved in the mixed model is significant, of the 
order of 100, and follows quantitatively well the open clus- 
ter observations, such as in the Hyades, or older clusters 
such as NGC 752. On the other hand, model ref burns its 
lithium mainly in the PMS phase with only a tiny fraction 
being depleted during the main sequence evolution due to 
gravitational settling. 

In Fig. [l] we also compare the sound speed and den- 
sity profiles of a typical static model which has been con- 
structed using the same physical inputs as in model B tz 
including its composition profiles as well. In the interior, 
the sound speed of the static model is almost the same as 
that of the evolutionary model, while in the outer layers 
it represents better that of the Sun. This improvement is 
most probably due to better treatment of surface layers in 
the static model through the use of a different prescrip- 
tion to calculate convective flux and also the adoption of 
a better atmospheric model. 

Despite all the positive aspects, solar models including 
tachocline mixing still need further improvements. Model 
B tZl for example, shows significant departure from the 
seismic Sun in its density profile, even though this quantity 
is in better agreement in comparison to model ref. Further, 
it assumes a photospheric value for the heavy elements, 
(Z/X) s ~ 0.0255 which is somewhat higher than the ob- 
served value. The reason is that a model calibrated to 
{Z/X) s = 0.0245, such as model B t in Brun et al. (1999), 
is not as close to the seismic Sun as model B tz , because of 
a smaller Z content and to the resulting variation of the 
opacities in the radiative interior. Therefore, we would like 
to reach between our new models and the Sun an agree- 
ment in density and sound speed better than that with 
model B tz , but without having to relax the (Z/X) s con- 
straint in the calibration process. Finally, the lithium de- 
pletion in the PMS phase is overestimated, indicating the 
need for a better treatment of this early phase of evolu- 
tion. All the cited improvements can come both from a 
better treatment of the mixing at the top of the radia- 
tion zone or by a better microscopic description. We refer 
to Piau & Turck-Chieze (2001) for a careful study of the 



PMS lithium depletion problem and focus our attention 
on the present Sun. We propose to keep for most of our 
models the same treatment for the tachocline mixing as in- 
troduced in Brun et al. (1999) but to allow for variations 
within uncertainties of the main physical ingredients, in 
order to see if any improvements can be obtained before 
introducing a new description of the tachocline mixing. 

3. Improved solar models 

3.1. Model parameters 

We see from the foregoing discussion that there is still need 
to improve our mixed tachocline models. In this section we 
outline the modifications and models computed for this 
study, which for the sake of clarity are also summarised 
in Table with their designation and the corresponding 
choice of parameters. 
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Fig. 2. Contribution to the nuclear energy generation of 
the main pp-chain and CNO cycles nuclear reactions as a 
function of the normalised solar radius. 

• pp nuclear reaction: When dealing with the prop- 
erties of the solar core one is inevitably led to the nu- 
clear reaction rates and their intrinsic uncertainties (Brun 
et al. 1998, Bahcall et al. 1998a; Morel et al. 1999). 
Unfortunately, the fundamental nuclear reaction in the 
pp-chain, i.e., p + p — > D +e + + v, is an clectroweak in- 
teraction and its cross section has not been determined 
via direct experimentation. Recent theoretical works give 
an uncertainty of about 1-2% in the determination of Sn 
(Adelberger et al. 1998). But a seismic calibration of this 
cross section indicated that up to a 4% increase is favoured 
(Antia & Chitre 1998). We have therefore decided to allow 
in our models for a variation up to 4% of this dominant 
cross-section and will also determine which value of S'ii(O) 
gives the best agreement with helioseismic inversion. 

Aside the fundamental pp nuclear reaction, other re- 
actions can possibly modify the central structure of the 
Sun. In Fig. | we plot the energy production of the dom- 
inant nuclear reactions from the pp chains as well as 
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Note: The parameters, d and 7V/27T, represent twice the extent of the tachocline h and the Brunt- Vaisala frequency (in fj,Hz) 
in the overshoot region and are related to our effective macroscopic coefficient Dt (cf., Brun et al. 1999). (Z/X) a , S11, S33, 
S34, fsc and Di are respectively, the surface ratio of heavy elements to hydrogen abundances of the models at the solar age, 
the variation applied in the nuclear cross section of pp, 3 He— 3 He and 3 He— 4 He, the screening prescription used in the model 
(either Weak or Intermediate) and the variation of the microscopic diffusive coefficient. 



the three main reactions from the CNO cycle. These re- 
actions are the following: D(p, 7) 3 He, 3 He( 3 f/e, 2p) 4 He, 
3 He(a,7) 7 Be, 7 Li(p, a) 4 He from the pp chains and 
13 C(p, 7 ) 14 N, 14 N(p, 7 ) 15 0(e+^) 15 N, 15 N(p,a) 12 C from 
the CNO cycle. To plot these curves we have used the 
internal structure of ref model for the present Sun. 

Not all the nuclear reactions displayed in Fig. || are 
expected to modify the thermal structure of the solar core. 
We can already discard some of them by considering either 
their importance in the energy budget or the time they 
require to reach equilibrium. 

• CNO cycle nuclear reactions: Because the Sun is a 
low mass star, we don't expect the CNO nuclear reactions 
to significantly influence the solar central region since the 
cycle contributes to less than 2% of the total nuclear en- 
ergy production (Clayton 1968; Bahcall 1989; Bahcall, 
Pinsonneault & Basu 2001). We have therefore not intro- 
duced any modification of these reactions in our models 
and refer to Turck-Chieze et al. (2001) for a discussion of 
their influence on the solar structure. 

• D-p and 7 Li-p nuclear reactions: Because the life- 
time of deuterium and lithium is very short, these ele- 
ments quickly reach equilibrium abundance in the tem- 
perature and density range prevailing in the solar core. 
Consequently, these reactions do not have a significant 
impact on the core structure even though there are very 
energetic. We refer to Gautier & Morel (1997) for a dis- 
cussion of the important D/H astrophysical ratio and to 
Brun et al. (1999) for a study of the influence of the 7 Li-p 
cross section on photospheric lithium abundance. We will 
adopt in all our models the cross sections proposed by the 
NACRE compilation (Angulo et al. 1999). 

• 3 He— 3 He and 3 He— 4 He nuclear reactions: These nu- 
clear cross sections play a crucial role in determining the 
branching between the ppl and the ppll and pplll chains 



and thus directly influence the high energy neutrino pro- 
duction. 3 He— 3 He is one of the most energetic reactions 
in the pp-chain along with T)-p. Unlike deuterium and 
lithium, 3 He does not reach its equilibrium value on a 
very short time scale but instead slowly builds up for 
temperatures less than 8 xlO 6 K (Clayton 1968). In a so- 
lar model, the resulting theoretical 3 He abundance profile 
peaks around r — 0.28i?o with a characteristic bell-like 
curve due to the competition respectively between its cre- 
ation and its destruction in the outer and in the inner 
regions of the solar core. At the same time, being one of 
the most energetic reactions in the pp chains and reach- 
ing equilibrium gradually in the outer parts of the nuclear 
region, this cross section is expected to have an influence 
on the thermal structure of the solar core at a level where 
seismic inversions can detect it. Even though the 3 He— 4 He 
nuclear reaction does not contribute much to the solar en- 
ergy budget, the fact that it involves both 3 He and 4 He 
chemical elements makes it also an important reaction to 
study. The experimental uncertainty of S33 and S34 are 
respectively ±8% and ±10% (Adelberger et al. 1998). 

• Intermediate screening: We use the intermediate 
screening prescription of Mitler (1977) in all our mod- 
els except one that uses the classical weak screening of 
Salpeter (1954). It should be recognised that the screening 
in stellar nuclear reaction rates is a sensitive issue which is 
not yet completely understood (Dzitko et al. 1995; Wilets 
et al. 2000). Depending on the solar thermodynamical con- 
ditions and the chemical species interacting in the nuclear 
reaction considered, one has to introduce the adequate 
screening factor, f sc , coming from the surrounding par- 
ticles present in the solar plasma which in general differ 
from the screening effect evaluated by nuclear physicist in 
their experiments. There are thus several sources of un- 
certainties in evaluating the cross sections and screening 
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effects for any given nuclear reaction that we intend to 
consider in this work. 

• Opacities, microscopic diffusion and heavy elements 
abundance Z: The structure of solar radiative zone is very 
sensitive to these two physical processes and to the heavy 
elements abundance. These are all closely related since a 
change in Z leads to a change in the opacity k and in the 
microscopic diffusion, which in turn modify Z as a con- 
sequence of the iterative calibration process. The opac- 
ity is accurately computed (error ~ 5%) for temperatures 
greater than 10 K (Iglesias & Rogers 1996), their main 
sources of uncertainties come from the relative composi- 
tion and the ionisation degree of the heavy elements as well 
as quantum effects (Rogers & Iglesias 1998). In this work 
we will leave k unchanged and will concentrate instead on 
the heavy elements abundance Z and on the amplitude 
of the microscopic coefficients Di. From a detailed com- 
parison with the work of Turcotte et al. (1998), Brun et 
al. (1998) have confirmed that the analytical expressions 
for the microscopic diffusion coefficients given by Michaud 
& Proffitt (1993) are accurate enough to deal with the 
solar case. For the Sun, these uncertainties are at most 
15%. In order to reduce the chemical composition gradient 
present at the base of the convection zone, which has been 
found to be too large in model ref compared to the seismi- 
cally deduced one, we have run one model with Di reduced 
by 10%. The heavy elements abundance of our models 
has been calibrated to the observed value of Grevesse et 
al. (1996), i.e., (Z/X) pb = 0.0245 (see Table g). Model 
B tz introduced earlier has been computed with an ini- 
tial heavy elements abundance Zq = 0.01959 = Z re f that 
leads to (Z/X) s = 0.0255 or 4% higher than (Z/X) ph . 

• Mixing: It is very tempting to introduce a mild mix- 
ing in the nuclear region, for example very close to the 
3 He peak, in order to improve the solar model structure 
and to reduce the predicted neutrino fluxes (Haxton 1997, 
Brun et al. 1998). However, the presence of a mixing in the 
solar core can be rejected on account of the helioseismic 
constraints, due to the huge disagreement in the central re- 
gion that it generates (of the order of few % in 5c/ c). This 
seismic evidence along with the recent results of the SNO 
neutrino experiment (Ahmad et al. 2001, 2002), strongly 
disfavour a "macroscopic mixing" as the source of the elec- 
tron neutrino deficit seen on Earth's detectors but instead 
support the idea of neutrino flavour oscillations. 

Therefore we limited ourselves to the tachocline re- 
gion for which we have better evidence for mixing and 
a relatively more elaborate physical description available. 
We used the parameters of model B tz (e.g., d = 0.10i?© 
and J\f/2n — 25 /iHz), which have been proven to give a 
reasonable agreement with seismic constraints and light 
elements photospheric abundance (cf. §2.3 and Brun et 
al. 1999). For one case we have assumed a wider mix- 
ing zone, i.e., d = 0.15Rq, in order to limit even more 
the gravitational settling and therefore reduce the steep 
gradients seen in the sound speed, density and composi- 
tion profile at the base of the convection zone. We have 
also used the prescription introduced by Elliott & Gough 



(1999), namely a constant diffusion coefficient Dte op- 
erating over a small domain d — 0.04i? Q . It is not our 
intention to reevaluate this coefficient, but just to com- 
pare both mixing prescriptions. 

Having introduced all these modifications in our evo- 
lution code we derive the profiles of sound speed, density, 
hydrogen abundance and temperature that we discuss in 
the following subsections. 

3.2. Inferred sound speed and density profiles 

We first consider the relative differences in the sound 
speed and density between the Sun and our new modi- 
fied solar models. 
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Fig. 3. The relative difference in sound speed and density 
profiles between the Sun and solar models as a function 
of fractional radius. Model ref including only microscopic 
diffusion and the mixed models B tz , NO and N (the lat- 
ter two with increased pp cross section), are represented 
respectively with solid, short-dashed, dash-dotted and long- 
dashed lines. Superimposed on model ref is the la error 
envelope coming from the helioseismic inversion. 



3.2.1. Increasing the pp cross section 

In Fig. U we display two new models, NO and N, which 
include respectively an enhanced pp nuclear cross section 
by 2% and 3.5%, along with the older models ref and B tz 
already presented in Fig. [l[ Clearly the relative differences 
in the sound speed and density between the Sun and model 
N, and to a lesser extent model NO, are smaller in com- 
parison to model ref and the mixed model B tz . This is 
really encouraging because contrary to model B tz , both 
models NO and N have been calibrated to the exact value 
of (Z/X) s — 0.0245 used in model ref. This indicates that 
the effect of Z via the opacity k can be compensated by 
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a small variation of the Sn. Thus a variation of the order 
of a few percent of either the heavy elements abundance 
or the fundamental pp nuclear reaction rate seems to have 
the same effect on the model sound speed profile in the up- 
per part of the radiative region but not below r — 0.3Rq 
where contrary to case B tZl model N departs from ref. 
For the density, a variation of Sn modifies the profile ev- 
erywhere resulting in a significantly better agreement for 
case N, at least above r = 0.2Rq. Thus it is very useful 
to assess the accuracy of a model by considering both the 
sound speed and the density profiles, because it allows one 
to distinguish the impact of different physical processes on 
solar structure. In this particular case an increase of Sn 
by 3.5% seems to be favoured by helioseismology, as op- 
posed to an increase of 5n by 2% or of (Z/X) s . 

It may be noticed that in Fig. 3 the agreement in sound 
speed inside the convection zone is much better as com- 
pared to that for the same models in Fig. 1. This improve- 
ment arises because we have scaled the solar radius in the 
models by a factor of 1.0003 before forming the difference 
with the seismically deduced sound speed and density pro- 
files. This scaling of the radius appears to remove most of 
the discrepancy in the upper convection zone, but does not 
affect the interior. This may be expected since a correc- 
tion of 0.03% in radial distance is only a small fraction of 
scale height in the interior, while it can become compara- 
ble to the scale height in the photospheric layers, resulting 
in significant differences in outer regions. A better agree- 
ment with the Sun is also obtained in the outer convective 
zone in the case of a static model (Fig. 1), which assumes 
a different treatment of the surface layers. It would appear 
that uncertainties in treatment of these layers are respon- 
sible for the discrepancy in the outer convection zone. In 
standard solar models the surface is normally defined as 
the layer where the temperature equals the effective tem- 
perature. Because of significant uncertainties in treatment 
of surface layers, the position of the surface may not be 
correctly estimated in a solar model. Thus, we believe that 
the scaling of radius effectively corrects for this error. In 
all subsequent figures we have used this scaled radius when 
comparing the solar models with profiles inferred from in- 
versions. 

3.2.2. Modifying the 3 He- 3 He and 3 He- 4 He cross 
sections 

We have just seen that a small change in Sn can signif- 
icantly improve the agreement between solar models and 
the seismic Sun when combined with a treatment of the 
solar tachocline. In the same spirit we have computed a 
sequence of models including modifications of S33 and S34 
nuclear reaction cross sections (cf. § |3.l| ). Models Nl and 
N2 share the same tachocline macroscopic treatment, cal- 
ibration of heavy elements to (Z/X) s — 0.0245 and in- 
crease of S11 by 3.5% as model N, but differ by having 
respectively an increase of 533 by 8% and a decrease of 
SW by 10% (see Table §. 



0.002 




Fig. 4. The relative difference in sound speed and density 
profiles between the Sun and solar models as a function 
of fractional radius, showing the effect of modifying the 
cross sections S33 and 6*34. The mixed models N, Nl and 
N2 are represented respectively with solid, short- dashed 
and long-dashed lines. Superimposed on model N is the 
la error envelope coming from the helioseismic inversion. 

In Fig. |] we plot the relative differences in density 
and sound speed between the Sun and models N, Nl and 
A^2. We first notice that the applied modifications of the 
cross section 6*33 and S34 improve the core structure both 
in density and sound speed. A variation of S34 by —10% 
seems to affect more the very central region than a varia- 
tion of S33 by +8% does, even if we take into account the 
fact that 5*33 has been varied by a smaller amount. For 
model N2 the agreement in 8c/ c in the core improves by 
a factor 2 with respect to model N, whereas for model Nl 
it does only by 20% or so. For Sp/p, the influence of these 
two cross sections is more modest and results in a small 
gradual change of the profile over the solar radius. 

We now turn to models A^02 and N03 displayed in Fig. 
H along with model NO as reference. Models A^02 and N03 
are identical to model NO except that they respectively 
include a decrease of S34 by 10% and the cumulative op- 
posite variations of the nuclear cross sections 633 by —8% 
and S34 by +10% (see Table|), in order to modify the pp 
branching ratio such as to increase the high energy neu- 
trino flux. Model A02 exhibits a better core profile than 
model NO does by having both relative differences closer 
to zero. However, these two models do not differ from each 
other as much as their counterpart models N and N2 do 
and are within the la error bar. In overall the decrease of 
S34 represents in this case a small progress toward a better 
agreement with the seismic data. On the contrary, model 
A03 that has been computed on purpose with opposite 
variations of S33 and S34 is almost everywhere in better 
agreement than model AT0 is. Its density profile is signif- 
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Fig. 5. The relative difference in sound speed and density 
profiles between the Sun and solar models as a function of 
fractional radius, showing the effect of modifying the cross 
sections S33 and 534. The mixed models NO, N02 and 
iV03 are represented respectively with solid, short-dash- 
long-dashed and dash-dotted lines. Superimposed on model 
iVO is the lfi error envelope coming from the helioseismic 
inversion. 



0.002 
0.001 



0.01 






' 1 ' 1 1 ' ' ' 1 ' ' 

N2 

N02 

NOW 'jSj&\ 






%/ \ 


V 


H 1 ~\ 1 1 1 1 1 1 1 1 1 1 l_v h, 





0.2 0.4 0.6 0.8 

r/R e 



Fig. 6. The relative difference in sound speed and density 
profiles between the Sun and solar models as a function 
of fractional radius, showing the effect of changing the 
nuclear screening. The mixed models N2, N02 and NOW 
are represented respectively with solid, short-dashed and 
long-dashed lines. Superimposed on model N2 is the la 
error envelope coming from the helioseismic inversion. 



icantly closer to the Sun. The main exception is in the 
very central part of the solar core, where the 5c/ c is quite 
off and the Sp/p exhibits a pronouncedly curved shape. 
Model iV03 is an interesting solar model but our variation 
of the S33 and 634 cross sections is certainly too large and 
goes in the wrong direction for the central parts. 

In summary, these new sets of results confirm that it 
is possible to improve the overall agreement between the 
models and the seismic Sun by modifying, within their un- 
certainties, the rates of important nuclear reactions such 
as S33 and S34. Here it appears that the seismic data 
favour a decrease of S34 and an increase of S33 rather 
than the opposite. 

3.2.3. Nuclear screening 

We would like now to characterise the influence of the 
nuclear screening on the solar model sound speed and 
density profiles. We have therefore changed the screening 
prescription from intermediate to weak (cf. §3.1), on the 
top of all the modifications already introduced in one 
model, that we have chosen to be 7V02, and called that 
new model NOW. In Fig. ^ we compare results obtained 
for models N2, N02 and NOW, in order to quantify the 
influence of the nuclear screening on the solar core. In the 
very central part of model NOW the agreement seems 
to not be as good as in model N02. The relative sound 
speed difference profile between the Sun and model NOW 
for r < 0.2Rq reaches a value of the same order as model 



NO in Fig. ^| and thus the change of nuclear screening 
seems to compensate the change by 10% in S34 made 
in model N02. The density profile of model NOW is in 
better agreement with the Sun than model N02 above 0.3 
Rq. We should be careful in evaluating the effect on the 
core structure of using a different screening prescription, 
because model NOW includes an increase of Sii(0) 
already calibrated on the intermediate screening, but the 
inversion seems to indicate that a smaller increase of Sn 
is favoured if one has to use a weak nuclear screening 
instead (cf., Table 3). 



3.2.4. Microscopic diffusion 

As already stated in §3.1, we are interested in a reduction 
of the steep composition gradient at the base of the con- 
vection zone, which implies a decrease of the microscopic 
diffusion coefficients Di. Model ND includes such a de- 
crease by 10% of Di along with other variations identical 
to model N2. It is obvious from Fig. that a reduction of 
microscopic diffusion is not appropriate, since model ND 
systematically departs more from the seismic Sun in the 
radiation zone than model N2. It is well known that mi- 
croscopic diffusion modifies and improves the stratification 
in the radiation zone (Christenscn-Dalsgaard et al. 1993). 
Based on the sound speed and density profiles, an increase 
of Di is favoured rather than a decrease as in model ND. 
However an increase of Di would make the composition 
gradient at the base of the convective zone even steeper. 
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Thus it seems quite unlikely that the microscopic diffusion 
is the remaining source of the discrepancies seen in a solar 
model like N2. 
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Fig. 7. The relative difference in sound speed and den- 
sity profiles between the Sun and solar models. The 
mixed models N2, NM, NE and ND are represented re- 
spectively with solid, short-dashed, long-dashed and dash- 
dotted lines. Superimposed on model N2 is the la error 
envelope coming from the helioseismic inversion. 



3.2.5. Adjusting the tachocline mixing 

In the sequence of models we have just discussed, we were 
only concerned with the influences on the solar structure 
of variations in the nuclear and atomic input data. We are 
now interested in assessing the impact on 5c/ c and 5p/p 
of varying the macroscopic parameters used for the solar 
tachocline. In Fig. (?] we also display models NM and NE 
which include the same microscopic ingredients as model 
N2 but with respectively, a broader mixing region with 
d ~ O.15i?0 (i.e., tachocline thickness h ~ 0.075i?©) and 
another macroscopic treatment for the tachocline (Elliott 
& Gough 1999), over a shorter distance h ~ O.O2i? (see 
Table ||). We clearly see that model NE does not reduce 
the bump in 5c/ c or the steep gradient in 5p/ p at the base 
of the convection zone seen in model ref in Fig. as much 
as the two other models N2 and NM. Even though model 
NE includes the same nuclear cross section modifications, 
it is further away from the seismic Sun. One reason for 
this poorer agreement is that the tachocline mixing is too 
shallow. The value of 0.02 i? Q adopted in model NE for 
the extent of the mixing has been calibrated by Elliott 
& Gough using a static model of the Sun at the present 
age and by convolving with the inversion kernels after- 
ward. This convolution procedure using the inversion ker- 
nels makes the effective thickness of the mixing broader, 



mainly because of their overlapping radial resolution, and 
as a consequence Elliott & Gough found that ~ 0.02i? Q 
was large enough to get rid of the bump in 5c/ c. It is 
quite puzzling that the introduction of their mixing in our 
evolutionary model does not give at all the same result. 
Indeed model NE is significantly different from models 
TV 2 or NM, that we believe are in quite good agreement 
with the seismic Sun. It can be seen in Fig. ^| that an in- 
crease of the extent of the mixing region in model NM 
does not affect that much the density and sound speed 
profiles except for the slight modification close to the base 
of the convection zone. We therefore conclude that not all 
prescriptions for 'tachocline mixing' give the same result 
and that ours is quite efficient in suitably modelling this 
transition region. 

3.3. Hydrogen abundance profile and photospheric 
composition 

We have so far addressed the question of the influence of 
mixing and microscopic variations on the primary inver- 
sion quantities such as the sound speed and the density. 
We would now like to assess what are the consequences 
of such changes on the profile and photospheric value of 
the hydrogen abundance X p h. Using the secondary inver- 
sion procedure introduced by Antia & Chitre (1998), we 
have compared the hydrogen abundance profile in the Sun 
with that predicted by our new set of solar models. By as- 
suming the Z profile of the models we have calculated the 
difference of X between each model and the Sun. 

In Fig. H we represent the absolute difference in 
hydrogen abundance profile 5X between the solar models 
ref, B tz , NO and N (cf., Table 2) and the Sun. As with 
5c/ c and 5pj p displayed in Fig. ^, it is quite clear that the 
mixed models are closer to the inferred solar hydrogen 
abundance than the purely microscopic model ref is. 
This improvement occurs mainly close to the base of the 
convection zone and in the convection zone itself (i.e., 
indicating a closer photospheric value) as expected by 
the introduction of our shallow tachocline mixing. The 
reason for such an improvement is twofold: 

• Firstly, by introducing a macroscopic mixing at the 
base of the convection zone, we hinder the gravitational 
settling of the chemical species and as a result there 
are relatively more helium and heavy elements in the 
convection zone, thus reducing the hydrogen contribution 
in the plasma composition mixture to a value closer to 
the seismically inferred one. 

• Secondly, the existence of an extended/mixed plateau 
of the chemical composition due to the presence of a 
macroscopic mixing at the base of the convection zone is 
in better agreement with the seismically inferred X profile. 

This result, along with the improved sound speed and 
density profiles discussed in the previous subsection, con- 
firms the presence of macroscopic mixing at the base of the 
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Fig. 8. The absolute difference in hydrogen abundance 
profile between the models and the Sun is shown as a 
function of fractional radius. Model ref including only mi- 
croscopic diffusion and the mixed models B tz NO and N 
are represented respectively with solid, short- dashed, dash- 
dotted and long-dashed lines. Superimposed on model ref 
is the la error envelope coming from the helioscismic in- 
version. 

convection zone and the necessity to introduce this process 
in solar models. Figure || also reveals that a small increase 
by less than 4% of the cross-section of the fundamental 
nuclear reaction pp is favoured as well. Indeed models NO 
and N are significantly closer to the inferred hydrogen 
abundance over the whole radiative interior (except for a 
small region between 0.35 and 0.5 i?©) than either of the 
two other models shown. Further, Fig. || shows that the 
photospheric hydrogen content of model N (AT p h — 0.7333 
see Table 3) and model iVO are closer to the seismically 
inferred value of 0.732 ± 0.001 than both models B tz and 
ref. Thus the effect of varying the nuclear cross section Sn 
is to change the hydrogen surface abundance in the model 
via the calibration procedure, the maximum amplitude of 
such modification being located in the central region, lead- 
ing to an improvement of the agreement between the Sun 
and model N by at least a factor 2. 

While determining the X profile through inversions 
we also get an estimate of pp reaction cross-section, Sn 
that is required to match the observed solar luminosity. 
Table 3 lists the values obtained by assuming the Z profile 
of each of the models considered in this study. It is clear 
that this estimate is not sensitive to other properties of 
the model, except for the treatment of the plasma screen- 
ing. Beside model NOW, all other models yields a value 
Sn w 4.06 x 10~ 25 MeV Barns, which is 1.5% higher than 
the value given by Adelberger et al. (1998). This is some- 
what less than the value given by Antia & Chitre (1999) 
or by Degl'Innocenti et al. (1998). This difference is due 



Table 3. Seismic inference on hydrogen abundance and 
S n 



Models 


X V, 


X- 




Rp f 


fl 73Q9 


n 7^1 1 




B tz 


0.7304 


0.7269 


4.053 


NO 


0.7338 


0.7322 


4.066 


N02 


0.7339 


0.7322 


4.054 


N03 


0.7334 


0.7323 


4.085 


NOW 


0.7337 


0.7322 


4.017 


N 


0.7333 


0.7323 


4.067 


Nl 


0.7334 


0.7323 


4.060 


N2 


0.7336 


0.7323 


4.054 


ND 


0.7321 


0.7327 


4.057 


NM 


0.7329 


0.7324 


4.055 


NE 


0.7362 


0.7317 


4.048 



Note: X p h and Xj nv correspond respectively to the photo- 
spheric hydrogen abundance achieved in the model and de- 
duced by seismic inversion. The seismically deduced pp reac- 
tion cross-section Sn for each models is given in unit of 10 -25 
MeV barns. 

to different treatment of the plasma screening in calcu- 
lating nuclear energy generation rates. If one uses weak, 
intermediate or strong screening for the solar plasma, the 
resulting increase in the pp cross section Sn found by 
seismic inversion will vary respectively between 0.5% and 
4%. Anyway all screening treatments seem to indicate a 
higher value of Sn than currently calculated by nuclear 
physicists. However, it appears that with our intermediate 
screening (e.g., Mitler 1977), an increase of Sn by 3.5% 
pushes the value of the pp cross section beyond what is re- 
quired from seismic constraints, but nevertheless we keep 
these models as they amplify the effect of increasing Sn 
and are in rather good agreement with the primary inver- 
sions of the sound speed and density. 

Apart from Sn Table 3 also gives the seismically in- 
ferred photospheric hydrogen abundance (X- mv ). Again 
this value is not sensitive to small differences in the Z 
profile, but is mainly determined by the photospheric Z 
value. Thus all models other than B tz g ive X[ nv ~ 0.732 
which is only slightly less than the value in the corre- 
sponding solar model (Xph). Further, this value will yield 
helium abundance, Y ~ 0.25, which is also close to the in- 
dependently inferred value using seismic inversions in the 
convection zone (Basu 1998, DiMauro et al. 2002). Thus 
the photospheric helium abundances obtained using dif- 
ferent techniques are consistent with each other. 

We now consider the effect of a variation of the 3 He- 
3 He and 3 He- 4 He reaction rates on the hydrogen abun- 
dance via their influence on the creation and destruction 
of helium in the solar core. In Fig. ^ we display the abso- 
lute difference SX between the mixed models N, Nl and 
N2 and the Sun. As we previously did, we prefer to use the 
mixed model N for comparison in this plot, since model N 
is significantly closer to the seismic Sun than model ref is. 
The absolute differences SX for the three models shown 
are all very close to zero. We find that as with the sound 
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Fig. 9. The absolute difference in hydrogen abundance 
profile between solar models and the Sun as a function 
of fractional radius. The mixed models N , Nl and N2 
are represented respectively with solid, short- dashed and 
long-dashed lines. Superimposed on model N is the la er- 
ror envelope from the hclioscismic inversion. 
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Fig. 10. The absolute difference in hydrogen abundance 
profile between solar models and the Sun. The mixed mod- 
els AT), N02, N03 and NOW are represented respectively 
with solid, short-dashed, long-dashed and dot-dashed lines. 
Superimposed on model NO is the la error envelope from 
the hclioscismic inversion. 



speed profile, variations of the nuclear reaction cross sec- 
tion S34 modifies relatively more the hydrogen abundance 
profile than variations of 633 , but obviously less than mod- 
ifications that variations of Su can produce. But for the 
quantity SX, model N gives a better agreement with the 
Sun than models Nl and N2. 

In Fig. [To] we display the difference in hydrogen abun- 
dance between the sequence of mixed models NO, N02, 
N03 and NOW and the Sun. A first small difference with 
the previous figure is that an increase of pp by +2% does 
not reduce as much the disagreement with the inferred hy- 
drogen abundance profiles than an increase by 3.5% does. 
As for model N2, the hydrogen abundance profile of model 
iV02 is only slightly affected in the inner central part. This 
seems to indicate that the variation of —10% of S34 is too 
large. Model N03, with opposite variations of £33 and S34 
cross sections, is shifted downward by —0.0005 compared 
to model NO over most of the radiative zone and its core 
profile is way off. All these models confirm the feeling that 
the effects of the S33 and S34 cross section are mixed, in 
the sense that in some regions the agreement is improved 
while in other regions it becomes worse. Thus it is diffi- 
cult to conclude if the variation of these cross sections is 
justified basing one's argument on seismic inversion of the 
hydrogen abundance. 

For model NOW, computed with a weak nuclear 
screening instead of intermediate as in model A^02, the 
result is a very small downward shift of 5X compared to 
model N02, but no significant improvement otherwise. It 
is thus unlikely that the screening effect will correct near 



the solar core the remaining discrepancy seen in the hy- 
drogen abundance profiles obtained in our mixed models. 

Finally in Fig. 11 we display the difference in hydro- 
gen abundance between the solar models and the Sun for 
cases NM and NE that include different macroscopic 
parameters and case ND with reduced microscopic dif- 
fusion along with case N2. For model ND the effect of 
the microscopic diffusion is subtle to be appreciated be- 
cause it modifies non uniformly the hydrogen abundance. 
The model seems to possess the closest hydrogen surface 
abundance A p h = 0.7321 relative to the inferred value 
-Xinv = 0.7327, and the smallest composition gradient at 
the base of the convection zone, thus justifying the use of a 
smaller microscopic diffusion coefficient Di . But the reduc- 
tion seems overestimated because model ND is the only 
solar model to exhibit a smaller hydrogen photospheric 
abundance X p h than the seismically inferred one. Further, 
deeper down in the radiative interior the hydrogen abun- 
dance profile departs too much from the inferred profile. 
Therefore, as for the sound speed, a reduction of Di seems 
to be discarded by present hclioscismic inversions. 

The effect of a broader mixing on SX (i.e., model NM) 
is to reduce the composition gradient and to extend the 
mixed plateau, properties that seem to be in better agree- 
ment with the inferred X profile. Moreover, X p h is closer 
to Xi n v in models NM and ND than for example model 
N2. Deeper down the improvement is not as obvious and 
considering the fact that the extent of the tachocline mix- 
ing is certainly over estimated in this model, we can hardly 
conclude that it constitutes a better solution than, say, 
models N or N2. Model NE is clearly worse in X as well 
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Fig. 11. The absolute difference in hydrogen abundance 
profile between solar models and the Sun as a function of 
fractional radius. The mixed models N2, NE, NM and 
ND are represented respectively with solid, short-dashed, 
long-dashed, and dot-dashed lines. Superimposed on model 
N2 is the la error envelope coming from the helioscismic 
inversion. 

as in 5c I 'c and 8p/p compared to other mixed models. 
Whereas its profile deep in the radiative interior is in rea- 
sonable agreement with the seismically inferred one, while 
it is clearly not the case for the upper part (r > 0.6i? Q ). 
The resulting composition gradient is too steep and the 
hydrogen abundance in the convective envelope too high. 

To summarise our findings about the hydrogen abun- 
dance profile, it can be stated that 'tachocline mixing' 
is very likely to occur in the Sun, and that an increase 
of Sn is clearly favoured by current helioscismic data, 
with an amplitude of the order of a few %. Less obvi- 
ous are the effects of the two others dominant nuclear 
cross sections, i.e., S33 and S34, although they can lead 
to some improvement as well. On the other side, neither 
the microscopic diffusion nor the screening effect seem to 
cause significant changes. However, the screening prescrip- 
tion has been found to modify significantly the value of 
the cross section of the pp nuclear reaction inferred by 
seismic inversion. Our last two models NE and NM in- 
cluding variations of the macroscopic parameters are not 
favoured by our study, because with model NE, the mix- 
ing is too shallow and does not vary with time and with 
model NM the mixing is too broad. We thus find that 
the seismically inferred photospheric hydrogen abundance 
is X[ DV = 0.732 ± 0.001, but the value change quite a bit 
depending on the value of heavy elements abundance as- 
sumed. 

We would like now to briefly discuss the new photo- 
spheric composition obtained by our models at the so- 
lar age, with particular emphasis on light elements deple- 



tion. As already stressed, the mixed models exhibit a bet- 
ter overall chemical composition, say compared to models 
with only microscopic mixing such as ref. One chemical 
element, namely the lithium, is crucial to assess the effi- 
ciency and time dependence of macroscopic mixing. With 
the presence of an effective macroscopic mixing at the base 
of the convection zone, all the models presented in this 
study are expected to deplete a fair amount of 7 Li. All 
models, except NE, burn indeed a substantial quantity of 
7 Li both in the PMS and in the main sequence phases, thus 
reaching a photospheric abundance at the solar age, Li s , 
significantly smaller than the initial/meteoritic value, Li^, 
i.e., Lio/Lig ~ 130 — 180, in reasonable agreement with 
the observations of Grevesse et al. (1996), Li /Li s ~ 140. 
Model NE depletes a large amount of lithium as well 
(Li /Li s ~ 100), but most of it (~ 90%) in the PMS 
phase, which is not realistic when compared with open 
cluster observations. This comes about because the macro- 
scopic coefficient Dte used in this model, following Elliott 
& Gough treatment of the tachocline, does not include any 
time dependence. Brun et al. (1999) have demonstrated 
that a proper time dependence of Dt causes significant 
lithium burning along the main sequence as well. Another 
important constraint is provided by the beryllium abun- 
dance, which requires that the mixing must be shallow in 
order not to destroy this element by more than 10% (see 
Brun et al. 1999; Bell, Balachandran & Bautista 2001). 
With either of the two diffusion coefficients used in this 
study, Dt or Dte, we easily achieve this goal, i.e., all the 
models but one deplete 9 Be by less than 10%. Model NM, 
which has been computed with the broadest tachocline 
mixing (i.e., h — 0.075Rq), leads to an underabundance of 
9 Bc of 20%, thus confirming that in this model the mixing 
extends too deep inside the radiative zone. We interpret 
this result as an indication that the tachocline mixing can 
not be broader than 5% in solar radius. Thus using this up- 
per limit for the tachocline extent, we are quite confident 
about the efficiency of our macroscopic time dependent 
coefficient DT{r,t) to model the tachocline region and to 
lead not only to the proper photospheric composition at 
the solar age but also in the earlier phases. However, as 
already stated with model Bt z , with our prescription the 
lithium depletion is still too big in the PMS, even though 
it is reasonably distributed over the whole temporal evo- 
lution compared, say, to model NE. 

3.4. Temperature profile and neutrino production 

It has been known for more than thirty years that stan- 
dard solar models and neutrino experiments on Earth dis- 
agree on the amount of neutrinos produced in the ther- 
monuclear core of the Sun, the former predicting always a 
flux in excess (Bahcall 1989). A number of ingenious sug- 
gestions have been given to either explain the discrepancy 
from revised and 'non-standard' solar models or by in- 
voking neutrino flavour oscillations from the electron neu- 
trino v e , generated in the pp chains and CNO cycles, to 
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its siblings the muon and tau v T neutrinos or to the so- 
called 'sterile neutrino' (Haxton 1995; Bahcall, Krastev 
& Smirnov 1998b). Here we intend to use the seismic di- 
agnosis to constrain as much as possible the theoretical 
neutrino flux, which is very sensitive to the central tem- 
perature. 

Figure |l2] represents the relative temperature differ- 
ence between solar models ref, B tz , N, NO, N02 and 7V03 
and the Sun. The overall agreement of the six models is 
quite satisfactory, with model ref being the least accurate 
in the tachocline region as expected. 
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Fig. 12. The relative temperature difference between so- 
lar models ref, B tz , N, NO, N02 and N03 and the Sun. 
Note the somewhat larger uncertainties for this variable 
in the neutrino production region (r < 0.3i?o) compared 
for example to the sound speed ones. 

We can indeed notice that the introduction of a macro- 
scopic mixing in the tachocline improves the profile of all 
the models by at least a factor 2 above r = O.6i?0, over 
model ref, thus confirming the importance of taking into 
account the mixing present in this transition region. In 
the bulk of the radiative zone, some improvement comes 
from the increase of the pp cross section and results in 
a flattening of ST/T and a slightly closer agreement with 
the Sun. For the very central part, where the neutrinos are 
produced, the secondary seismic inversion of the temper- 
ature is less accurate with an error bar of ±2 x 10 -3 , and 
thus does not constrain the solar models as much. Models 
ref and B tz are surprisingly good there. Model N is quite 
satisfying except in this very central part. This can be 
interpreted as an excessive increase of the pp cross sec- 
tion. Models NO and N03 seem to be our best models in 
the range [0,0.7] Rq. While this could have been expected 
from model NO, based on the inversion of the sound speed 
and density profiles (§3.2), it was not so for model N03 
that has been computed with opposite variation of the nu- 



clear cross sections 6*33 and 6*34. Effectively, it seems that 
for the temperature profile, the decrease of S34 by 10%, 
as in model iV02, does not lead to any progress in the 
central part while an increase does. This conclusion is at 
odd with what we learnt from the sound speed inversion. 
This could mean that either the secondary inversion of the 
temperature is not as reliable as the sound speed (that is 
in part true but unsatisfactory) or that the temperature 
varies differently with modifications of the main physical 
ingredients than the sound speed does, due for example 
to a compensatory change of the central composition. So, 
we have to be cautious in our conclusions regarding the 
very central part since this is the region where the inver- 
sions are the least reliable. Nevertheless, we still consider 
that for the temperature profile models ref, B tz , NO, N03 
and more marginally, model N , all represent seismically 
acceptable solutions of the solar core. As a result, what 
one can expect to be the impact on the neutrino fluxes of 
such diverse temperature profiles? 

Table 4. Neutrino Fluxes at Earth 







Detector 




Models 


71 Ga 


37 C1 


Water 


Ref 


127.1 


7.04 


4.99 


Btz 


127.1 


7.04 


4.99 


NO 


123.7 


6.41 


4.48 


N03 


128.2 


7.08 


4.99 


NOW 


121.4 


6.10 


4.25 


N 


122.5 


6.18 


4.29 


Nl 


121.3 


6.0 


4.16 


N2 


119.3 


5.7 


3.93 


ND 


118.7 


5.59 


3.85 



Note: The Ga, CI and Water columns correspond respec- 
tively to the predicted solar neutrinos fluxes for the Gallium, 
Chlorine and Water experiments (cf., Table [l]). The gallium 
and chlorine neutrino fluxes are given in SNU whereas the wa- 
ter (SNO/SuperKamiokande) ones are in 10 6 cm -2 s _1 . 

To answer that question we have summarised in Table 
[| the neutrinos fluxes of the most significant models. By 
comparing for example models ref and N2, we find that 
for the latter the Gallium flux is reduced by 8 SNU, the 
Chlorine by 1.34 SNU and the 8 B by 10 6 cm~ 2 s" 1 down 
to 3.93 x 10 6 cm~ 2 s _1 . Such theoretical fluxes are still too 
high compared to the neutrino experiment on Earth (cf., 
Table 1), if no other modifications are introduced either in 
the model or in the quantum properties of the neutrinos. 
For example, the production of the 8 B neutrinos is directly 
sensitive to the p- 7 Be nuclear cross section. The value used 
in this study is 5x7(0) = 19.1^2 eV barns (Adelberger et 
al. 1998), and corresponds to an intermediate value com- 
pared to the recent estimate of Davids et al. (2001) (i.e., 
5x7(0) = 17.8±\i eV barns) or of Junghans et al. (2002) 
(i.e., 5x7(0) = 22.3 ± 0.7 ± 0.5 eV barns). Thus we can 
conservatively consider that the error bar in 5x7 nuclear 
cross section is at least of the order of ±10%. Such an un- 
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certainty, results in an increase or decrease of our 8 B flux 
by ±10% as well. 

By taking into account most of the uncertainties 
present in a solar model to calculate the neutrino fluxes, 
such as the nuclear reactions cross sections, the screening, 
the heavy elements abundance, the absorption cross 
sections for Gallium and Chlorine experiments, the 
amount of microscopic and macroscopic diffusion, etc. 
(see Turck-Chieze et al. 2001), we end up with the follow- 
ing fluxes (errors have been obtained by quadratic sum of 
the individual contributions assumed to be independent): 

71 Ga= 123.7 ±8.7 SNU, 

37 C1= 6.41 ± 0.86 SNU and 

Water= (4.48 ± 0.71) x 10 6 cm" 2 s -1 . 

These values are a bit lower than Bahcall et al. (2001) 
and Turck-Chieze et al. (2001), but remain within the 
la error range. However, they are still significantly larger 
than the observations, unless one invokes non standard 
neutrino properties. Such evidence for neutrino flavour os- 
cillation have been recently given by Ahmad et al. (2001) 
based on the careful study of the high energy neu- 
trino fluxes detected by the SNO detector and by the 
SuperKamiokande experiments (Fukudact al. 1998). Since 
then SNO has also measured the 8 B neutrino flux using 
the neutral current channel, which is equally sensitive to 
all neutrino flavours (Ahmad et al. 2002). The resulting 
8 B flux, is found to be (5.09 ±0.62) x 10 6 cm~ 2 s -1 , which 
is within the error bars of the current solar neutrino pre- 
diction. Our helioseismic study seems to favour more the 
lower range of the detection than the upper range. 

4. Discussion and conclusions 

In this work our goal was to assess the effect of uncertain- 
ties in nuclear reaction rates, atomic data and diffusion 
coefficients on solar models, which we have compared with 
the results of helioseismic inversions. We have concluded in 
§[| that a variation of the nuclear reaction rates can have a 
significant impact on the solar structure and that current 
modelling coupled with seismic data favours some change 
of the accepted central value of three important nuclear 
reaction rates, pp, 3 He- 3 He and 3 He- 4 He. While the pp re- 
action rates needs to be increased by about 1.5% over the 
currently accepted value from Adelberger et al. (1998), the 
constraints on other reaction rates are less clear. Even the 
increase of pp rate involved in this study is less than that 
inferred by Antia & Chitre (1999) and Degl'Innocenti et 
al. (1998) and the difference can be attributed to different 
treatment of nuclear energy generation in CESAM as com- 
pared to the version of Bahcall's energy routine used in 
earlier estimates. More precisely, most of the discrepancy 
comes from the different screening prescriptions used. In 
the previous study of Antia & Chitre (1999), the screen- 
ing formulation of Graboske et al. (1973) was assumed as 
opposed to Salpeter (1954) or Mitler (1977) in this new 
study. We refer to Dzitko et al. (1995) and Wilets et al. 



(2000) for a detailed comparison of several weak, interme- 
diate and strong, screening prescriptions. 

Our work confirms the result of Brun et al. (1999) 
that the implementation of macroscopic mixing in the 
tachocline improves the agreement between solar mod- 
els and seismic Sun (§3.2). In particular, direct compar- 
ison of hydrogen abundance profiles between our mod- 
els and the Sun as inferred from seismic inversions has 
demonstrated that the models with tachocline mixing are 
in much better agreement with the Sun in tachocline re- 
gion. Comparison of hydrogen abundance profile in solar 
models with inferred profiles show that there is still a dis- 
crepancy of about 0.003 (§3.3). The largest discrepancies 
occur in the region close to the tachocline and at about 
r = 0.2_Rq, where the sound speed also shows maximum 
discrepancy. Thus it appears that there is still some scope 
for improving the formulation for calculating the mixing 
in tachocline region. We were not able to achieve any sig- 
nificant improvement compared to the results of Brun et 
al. (1999) by adjusting the parameters of the tachocline 
model. With some modifications the agreement improves 
in the tachocline region, like in model NM, but it tends to 
worsen in other places. We have found as well that based 
on the 9 Bc photospheric depletion, the maximum extent 
of the mixing in the tachocline is 5% of solar radius. It 
is quite possible that a major part of the remaining dis- 
crepancies (about 0.1% in sound speed, 1% in density and 
0.003 in X) in our improved models may be due to uncer- 
tainties in input physics, like the opacities or equation of 
state. 

Comparing the surface hydrogen abundance in solar 
models and those obtained by inversions (Table 3) with 
the abundance inferred directly from seismic data (Basu 
1998), we again find a difference of 0.001. This suggests 
that discrepancies of this order are present in the in- 
put physics which is used in inferring these independent 
measurements of surface hydrogen abundance. The he- 
lium abundance drawn from the structure of the ionisa- 
tion zones is particularly sensitive to the equation of state 
(Basu & Antia 1995), while its abundance inferred from 
the solar models depends on the whole input physics used 
in constructing those. Similarly, the seismically inferred 
hydrogen abundance profile is sensitive to opacity as well 
as to the Z profile used in inversion. Thus the discrepancy 
between these three independent estimates of surface hy- 
drogen abundance probably gives an estimate of errors in 
input physics. The difference between solar models with 
tachocline mixing and seismically inferred profiles is of the 
same order and could be due also to remaining uncertain- 
ties in input physics. 

As seen in Fig. 1 the evolutionary solar models show 
a significant departure from seismic inferences inside the 
convection zone. We find that scaling the radius of these 
models by 1.0003 before taking the differences with in- 
verted profile removes most of the discrepancy in the con- 
vection zone. We believe that this is due to uncertainties in 
the treatment of surface layers. In the solar models the sur- 
face is defined by the layer where the temperature equals 
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effective temperature and it is quite possible that because 
of uncertainties in surface layers the location of this point 
has an error of about 200 km in our evolutionary models. 
We would like to point out that adjusting the solar radius 
in the model by 200 km does not remove this discrepancy 
since even in the new model the position of the surface 
has the same uncertainty. Thus this scaling of radius has 
no relation to uncertainty in solar radius itself. 

The calculated neutrino fluxes in solar models with 
tachocline mixing are found to be somewhat lower than 
those in standard solar model of Bahcall et al. (2001), but 
within the la error limits (§3.4). The main reason for re- 
duction is the increase in pp reaction rate, which reduces 
the central temperature required to generate the solar lu- 
minosity. Recently, using both the charged and neutral 
current channels, SNO has measured the total 8 B neutrino 
flux of 5.09 x 10 6 cm" 2 s" 1 (Ahmad et al. 2002). This is 
somewhat larger than the value we find in our models. In 
the modified models Nl, N2 the neutrino fluxes are much 
lower, because we have chosen to modify the nuclear reac- 
tion rates to reduce these fluxes. If instead as with model 
A^03 we had modified S33 and S34 in opposite directions 
the neutrino fluxes would have increased and it would be 
possible to get values close to that inferred by Ahmad et 
al. (2002). Also our increase of pp reaction rate by 3.5% 
in some of these models is probably an overestimate, since 
seismic models with same input physics appear to need 
an increase by only 1.5% to produce the required solar 
luminosity. When almost all the sources of uncertainties 
are taken into account in evaluating the theoretical neu- 
trino fluxes, our result is within la of the other published 
theoretical neutrino fluxes and the recent SNO compila- 
tion. It will be interesting in the near future, when the 
SNO collaboration results will have been integrated over 
a longer time to see if the helioseismic tool will be able 
to constrain even more efficiently the main nuclear cross 
sections, via a careful analysis of the neutrino spectrum, 
and also delineate the parameter space for mass-squared 
difference-mixing angle plane. 
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